Evolution of the Spin Gap Upon Doping a 2-Leg Ladder 



O 

o 
o 



> 

o 

in 
o 
o 
o 



c 

o 
o 



13 



D. Poilblanc a *, 0. Chiappa a and J. Riera b t, 
a Laboratoire de Physique Quantique & UMR-CNRS 5626, Universite Paul Sabatier, F-31062 Toulouse, France 
Instituto de Fisica Rosario, Consejo Nacional de Investigaciones Cientificas y Tecnicas, y Departamento de Fisica, 
Universidad Nacional de Rosario, Avenida Pellegrini 250, 2000-Rosario, Argentina 

S.R. White c * and D.J. Scalapino rf § 
c Department of Physics and Astronomy, University of California, Irvine, CA 92691 
d Department of Physics, University of California, Santa Barbara CA 93106 

(February 1, 2008) 

The evolution of the spin gap of a 2- leg ladder upon doping depends upon the nature of the lowest 
triplet excitations in a ladder with two holes. Here we study this evolution using various numerical 
techniques for a t-t'-J ladder as the next-near-neighbor hopping t' is varied. We find that depending 
on the value of t' , the spin gap can evolve continuously or discontinuously and the lowest triplet state 
can correspond to a magnon, a bound magnon-hole-pair, or two separate quasi-particles. Previous 
experimental results on the superconducting two-leg ladder Sri4_ a; Ca 2 :Cu24 04i are discussed. 

PACS: 71.27. +a, 75.50.Ee, 71.10.-w, 75.40.Mg 



Studies of strongly-correlated electrons confined to 
two-leg ladders and described by t-J and Hubbard mod- 
els have provided important insights into the high T c 
cuprate puzzle. These models are known to exhibit a 
gaped spin liquid state at half-filling and upon doping 
to evolve into a Luther-Emery state characterized by 
d a; 2_j / 2-like pairing and 4fcp CDW correlations 
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key feature of the Luther-Emery state is the existence of 
a gap Ag in the excitation energy of the spin degrees of 
freedom. At half-filling, the spin gap Ag is set by the 
K = (7r,7r) magnon excitation energy Am which is of 
order J/2 for an isotropic Heisenberg ladder with a near 
neighbor exchange interaction J. However, as discussed 
by Tsunetsugu et. al. ||, there can be a discontinuous 
evolution of the spin gap upon doping. In particular, 
they note that a pair can be dissociated into two charge 
|e| and spin S =1/2 quasi-particles, and the low-energy 
continuum for such scattering is set by the pair-binding 
energy Ap. Then, if the pair-binding energy is less than 
the half- filled spin gap (Ap < Am), there will be a dis- 
continuous decrease in the spin gap upon doping to a 
value equal to the pair-binding energy Ap. Thus, while 
there is still an S = 1, K = (tt,tt) magnon excitation 
with energy Am in the infinitesimally doped ladder, if 
Am > Ap a lower energy S = 1 state exists in which a 
pair is dissociated into two quasi-particles. 

There is a low energy continuum of excited states 
corresponding to two quasi-particles, each in an even 
parity k y = state, which have a total momentum 
(k x ,k y ) = (0,0). Here k y = for a bonding and it 
for an anti-bonding quasi-particle respectively. The sin- 
glet and triplet continua start at the same energy Ap. 
In addition to these scattering states, there can also be 
a bound S — 1 state in which a bonding and an anti- 
bonding quasi-particle with momentum k y — tt hybridize 
with a magnon excitation of the spin background [|],^| . 



If there is such a bound magnon-pair with energy Amp, 
then it will set the spin gap in the doped ladder provided 
Amp < Am |§- Such a scenario occurs e.g. in lad- 
ders with anisotropic rung (Jj_) and leg (J||) couplings 
in the range 0.4 < J±/J\\ < 1.4 0. Here we com- 
bine exact diagonalization (ED) Q and density-matrix- 
renormalization-group (DMRG) techniques Q to investi- 
gate the evolution of the spin gap when one pair of holes 
is added to a t-t'-J ladder. The next-near-neighbor one- 
electron hopping t' provides a useful tuning parameter 
to study the interplay of the magnon gap Am, the pair 
dissociation gap Ap, and the bound magnon-pair spin 
gap A mp in setting the spin gap Ag of the lightly doped 
ladder. 

The Hamiltonian for the t-t'-J ladder is 
H = J y^(Sj,> • Sj+i,;\ - ^n i;A n i+li A) (1) 

i 

+ t X^U^+bV + h - c -) + < X!( c !m c *,2, s + h.c.) 

i,\,s i,s 

Here cj , creates an electron of spin s on site i of leg 
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<C 4 ,A, S ')/ 2 and n i,\ 



Ci,\.s- We have taken both the near-neighbor 
leg and rung one-electron hopping matrix elements equal 
to t and the diagonal next-near-neighbor term equal to 
t' . The exchange interaction J is taken as isotropic be- 
tween near-neighbor leg and rung sites and throughout 
this J/t = 0.5. 

We begin with our conclusions shown in Fig. 1(a) 



1 



where we have plotted the excitation energies AE of var- 
ious triplet states versus t' . The spin gap A5 of the two 
leg t-t'-J ladder doped with two holes is defined as the 
difference between the ground state energies of the sys- 
tem with two holes and 5=1 and 5 = respectively. 

A s = E Q (n h = 2, 5 = 1) - E (n h = 2,5 = 0) . (2) 

The stars in Fig. 1(a) show As versus t'/t obtained from 
DMRG results on 2 x L ladders with L = 32. The dashed 
line is the DMRG result for the magnon excitation of the 
undoped ladder obtained from 

A M = E (n h = 0, 5 = 1) - E (n h = 0, 5 = 0) . (3) 

That this difference in ground state energies corresponds 
to the (71", 7r) magnon is known from ED calculations in 
which the momentum of the excitation is specified. The 
open diamonds show the triplet excitation energy in the 
K = (71", 7t) sector, obtained from a finite size scaling 
analysis using ED. Finally, the solid curve in Fig. 1 corre- 
sponds to the pair-binding energy calculated with DMRG 
from 
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Eo (n h = 2,5 = 0) +E (n h = 0, 5 = 0) 



2E ( n h = 1,5= - 



(4) 



with rih the number of holes relative to the half-filled 
ladder (in agreement with the ED results for t' = in 
Ref. 0). As shown in Fig. 1(b), Ap sets the two quasi- 
particle continuum. Here infinite size extrapolated ED 
results for the lowest energy excited singlet and triplet 
states in the K = (0, 0) sector are plotted as open sym- 
bols and the solid circles are DMRG data for the pair- 
binding energy Ap, Eq. (I|). These energies are in good 
agreement, consistent with a picture in which a pair dis- 
sociates into two quasi-particles. 

As discussed below, we have used ED, in which the mo- 
mentum of the state can be specified, as well as DMRG 
calculations of the hole and spin correlations in order to 
interpret the results shown in Fig. 1(a). Here we sum- 
marize what these show. Basically, there are three dif- 
ferent regimes set by t'/t. For —0.5 < t'/t < —0.2, the 
discontinuous drop in the spin gap with doping reflects 
the fact that the pair binding energy Ap is less than 
the (it, tt) magnon energy Am of the undoped ladder. 
Thus, when the system is doped, a singlet pair can dis- 
sociate into two separate quasi-particles with total spin 
5 = 1, reducing the spin gap A5 from Am to Ap. In 
this region, there is a bound magnon-hole pair with a 
minimum energy at (7r,7r) but its energy Amp is larger 
than Ap so that the spin gap is set by Ap. In the re- 
gion —0.2 < t'/t < 0.35, the situation changes. The 
pair binding energy Ap becomes greater than the energy 
to create a bound magnon-hole pair Amp, but Amp is 
less than the energy to create a separate magnon Am- 
Thus, in this parameter region the lowest energy triplet 
state of the 2-hole doped ladder has momentum (tt, tt) 



and corresponds to a bound magnon-hole pair so that 
A s = A M p- Finally, for 0.35 < t'/t < 0.5, the energy 
of the triplet K = (n, tt) excitation becomes equal to 
the 5=1 magnon energy of the undoped ladder. Here 
DMRG calculations of the spin and charge correlations 
show that the excitation corresponds to a magnon which 
is uncorrelated with the bound singlet pair. Thus, in this 
region, there is no discontinuity in the spin gap upon dop- 
ing. 
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FIG. 1. (a): Low energy triplet excitations of a 2-leg t-J 
ladder versus t' for J/t = 0.5. The shaded region corre- 
sponds to the onset of the two quasi-particle continuum set 
by the pair-breaking energy Ap. The dashed line denotes 
the magnon energy Am of the undoped ladder. The open 
diamonds show ED results for the lowest triplet excitation in 
the K = (tt, tt) sector and the stars denote the DMRG spin 
gap As measured for a ladder with 2 holes, (b): the open 
symbols show ED results for the lowest energy singlet and 
triplet excitations in the K = (0, 0) sector. The solid circles 
show DMRG data for the pair- binding energy Ap, Eq. (^J). 
The solid line through the DMRG points denotes the two 
quasi-particle continuum as also shown in (a). 

ED calculations were carried out on 2 x L ladders 
with L an odd number of sites. Both periodic and anti- 
periodic boundary conditions for L up to 13 were used 
fllPf . In Fig. 2 we show results for the triplet excita- 
tion energies in the K = (tt, tt) sector for a sequence of 
2 x L ladders with 2 holes for various values of t'. Here 
the excitation energy is measured relative to the 2-hole 
K = (0, 0) ground state. The lowest triplet K = (71", tt) 
state is found to be separated from a quasi-continuum 
of higher energy states. In Fig. 2, the error bars mark 
the difference between the results obtained using peri- 
odic and anti-periodic boundary conditions with the open 
symbols marking the mean value. Since the actual lon- 
gitudinal momentum for a finite ladder is 7r(l — 1/L), 
we have extrapolated these results using a scaling form 
A + B/L + c/L 2 . The solid symbols denote the DMRG 
calculation of the spin gap Ag, Eq. (||), for an open 
2 x 32 ladder with 2 holes (larger lattices are also in- 
cluded for t' = 0). For -0.2 < t' < 0.5, the extrapo- 
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lated ED results for the K = (%, tt) triplet pass through 
the DMRG spin gap Ag. However, for t' = —0.5, the 
DMRG determined spin gap lays well below the extrap- 
olated K = (tt,tt) triplet. As discussed in the introduc- 
tion, for —0.2 < t', the spin gap is set by the excitation 
in the triplet K = (tt, tt) sector. However, for t' < —0.2, 
the spin gap is set by the onset of the two quasi-particle 
continuum Ap which goes to zero as t' approaches —0.5. 
Note that, even when As < Amp, the magnon-hole pair 
state could still be locally stable if the decay process into 2 
quasi-particles with the same momentum (tt, 7r) is impos- 
sible. Although the ED results approach the K = (tt, tt) 
magnon energy of the undoped ladder for t' > 0.35, 
DMRG results show that the character of the triplet exci- 
tation changes from a bound magnon-hole-pair to a sep- 
arate magnon and hole-pair state. Thus, in this regime, 
the spin gap is set by the excitation energy of the magnon 
Am and is therefore continuous upon doping. 
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sponds to having the holes on diagonal sites. Note that 
the triplet-bound state is more extended than the singlet 
bound state. The lowest diagram shows the spin distri- 
bution for the S z = 1 triplet state when the two holes 
are projected out at their most probable sites. It is clear 
that this state is a bound magnon-hole pair. 
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FIG. 3. Structure of the ground state and the lowest energy 
triplet state of a 2-leg ladder with 2 holes and t' = 0. In the 
top and middle diagrams, the diameter of the black circles 
is proportional to the probability of finding the second hole 
when one hole is projected out at the center of the upper 
leg for the ground singlet and lowest energy triplet states 
respectively. Note that the triplet state is more extended 
than the singlet bound state. The bottom diagram provides 
the spin distribution in the triplet state. 



FIG. 2. A finite-size scaling analysis of the lowest triplet 
eigenenergies of a 2 x L ladder with 2 holes in the K = (tt, tt) 
sector measured with respect to the K = (0, 0) ground state 
singlet for various values of t' . The open symbols are ED data 
obtained on closed ladders up to 2 x 13 in size. The mean value 
between periodic and antiperiodic boundary conditions is de- 
noted by the open symbols. The error bars are given by the 
energy difference between the two boundary conditions. The 
full symbols are DMRG results for the spin gap As defined 
by Eq. (§). 

In order to get a clearer picture of the nature of the 
triplet excitations which determine the spin gap, we have 
used DMRG results to study the spin and hole correla- 
tions in these states. In Fig. 3 a center section of a 2 x 32 
ladder with t' = is shown. The upper and middle di- 
agrams show the probability of finding the second hole 
when the first hole is projected out at the center of the 
upper leg for the singlet (llj and triplet state respec- 
tively. In both of these states, the two holes are bound 
and the most probable configuration for J/t — 0.5 corre- 



Similar calculations for t' = —0.5 show that in the 
triplet state the two holes are unbound while for t' = 0.5 
the two holes are bound into a singlet and uncorrelated 
with the spin 1 excitation. This behavior is shown in 
Fig. 4, where we have plotted 

(5,(4)) = (S z (£ x , l)P h (i)P h (j)) I (P h (i)P h (j)) (5) 

versus £ x for t' — 0, —0.5, and 0.5. Here Ph(i) is the 
projection operator for a hole at the i th site. For t' = 0, 
we have set i = (16,2) and j = (17,1), corresponding 
to the most probable hole location. Here, as previously 
illustrated in Fig. 3, we see that the spin is bound to 
the hole-pair. For t' = 0.5 we again have a situation 
where the holes are most likely to sit close to each other, 
and here we have projected them onto i = (16, 1) and 
j = (16,2). However, in this case, the spin 1 is spread 
out corresponding to a magnon which is not bound to the 
hole-pair. Finally, for t' = —0.5, one finds that the low- 
est energy triplet excitation corresponds to two separate 
quasi-particles. 
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We finish with a brief discussion of some experimen- 
tal results for the superconducting two-leg doped ladder 
Sr 2 Cai2Cu2404i ]T3). Nuclear magnetic resonance mea- 
surements of the copper-63 Knight shift and relaxation 
time T\ suggest a collapse of the spin gap with pres- 
sure. We believe this signals the appearance of new low 
lying triplet excitations upon doping the ladder planes 
and points towards a negative value of t' Jll|. In this 
regime, due to the presence of a low-energy quasi-particle 
continum located predominantly around the zone center, 
momentum-resolved experiments like inelastic neutron 
scattering would be essential to search for sharp finite 
energy triplet excitations. 
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FIG. 4. Here, for three different values of t' we show 
{S x (£ x )} (as defined by Eq. (||)) versus l x . For t' = 0, we take 
i — (16,2), and j = (17, 1), one of the most probable loca- 
tions for the two holes, and one sees that a magnon is tightly 
bound to the hole pair. For t' = 0.5, we have i = (16, 1) 
and j = (16,2), also one of the most probable locations for 
the holes. Here we see a magnon which is not bound to the 
hole pair. For t' = —0.5, in order to illustrate how separate 
S — 1/2 spins are bound to each hole, we plot {S z (£x)) with 

1 = (8,1) and j = (25,1). 

To summarize, using ED and DMRG calculations, we 
have found that the spin gap can evolve in different ways 
when two holes are doped into a 2-leg ladder. When the 

2 holes are added, it is possible that the lowest energy 
triplet state simply remains the K = (tt, n) magnon so 
that there is no change in Ag. In this case the 2 added 
holes remain in a bound d x 2_ y 2-\ike singlet state and a 
triplet magnon similar to that of an undoped ladder is 
created. As the length of the ladder increases, the in- 
teraction between these two entities becomes negligible. 
We see this happening for i' > 0.35. It is also possible 
that the lowest energy triplet state has K = (0, 0) and 
is set by the two-quasi-particle continuum corresponding 
to the pair-binding energy Ap. In this case, there is a 



discontinuous change in the spin gap upon doping and 
the lowest energy triplet state arises from the dissocia- 
tion of a pair into two quasi-particles. We see this for 
the present model when t' < —0.2. Finally, for the inter- 
mediate region —0.2 < t' < 0.35 we find that the lowest 
energy triplet state has K = (71", tt) and corresponds to 
a bound magnon- hole pair with energy Amp < Am- In 
this case, there is again a discontinuous evolution in the 
spin gap from Am to Amp upon doping. 
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